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BIVARIATE CHEBYSHEV-I WEIGHTED ORTHOGONAL 
POLYNOMIALS ON SIMPLICIAL DOMAINS 

MOHAMMAD A. ALQUDAH 


Abstract. We construct a simple closed-form representation of degree-ordered 
system of bivariate Chebyshev-I orthogonal polynomials 9n,r (u, v, w) on sim- 
plicial domains. We show that these polynomials 9n,r(u , v,w), r = 0,1.. .. , n; 
n > 0 form an orthogonal system with respect to the Chebyshev-I weight 
function. 


1. Introduction and Motivation 

Orthogonal polynomials have been studied thoroughly mm, the Chebyshev 
polynomials of the first kind (Chebyshev-I) are among these polynomials. Although 
the main definitions and basic properties were considered many years ago pa ms, 
the cases of bivariate or more variables have been studied by few scholars. 

Chebyshev-I polynomials (T njr (u, v, w)) are defined in many textbooks U [8], 
they are orthogonal to each polynomial of degree < n — 1, with respect to the 
defined weight function. However, for r ^ s, T n ^ r (u,v,w) and T niS (u,v,w) are not 
orthogonal with respect to the weight function. 

Farouki [5] defined orthogonal polynomials with respect to the weight function 
W(u, B,«/) = lona triangular domain. These polynomials P n , r {u, v, w) are orthog¬ 
onal to each polynomial of degree < n — 1 , and also orthogonal to each polynomial 
P n ^ s (u,v,w), where r ^ s. 

In this paper, we construct bivariate Chebyshev-I weighted orthogonal polyno¬ 
mials ^ 7 n t r(u,v,w) with respect to the weight function W(u,v,w) = « _ 5ir5(l — 
w) J ,'y > 1, on triangular domain T. We show that these bivariate polynomials 
$n,r{u, v, w), r = 0,1 ,..., n; n > 0, form an orthogonal system over T with respect 
to W(u,v,w) = u~^v~ 2(1 — u>) 7 ,7 > 1. 


1.1. Bernstein and Orthogonal polynomials. Consider a triangle T defined 
by its three vertices p fe = ( Xk,yk), k = 1,2,3. For each point p located inside the 
triangle, there is a sequence of three numbers u,v,uj > 0 such that p can be written 
uniquely as a convex combination of the three vertices, p = upj^ + vp 2 + RP3, 
where u + v + w = 1. The three numbers u = areaf P.P2>P3) v = aT, ea(p 1 .p,p 3 ) 

area(p 1 ,p 2 ,p 3 ) ’ area{p 1 ,p 2 ,p 3 ) > 


2000 Mathematics Subject Classification. 42C05, 33C45, 33C70. 

Key words and phrases. Chebyshev, Orthogonal, Triangular, Bernstein, Bivariate Polynomials. 
©2010 Ilirias Publications, Prishtine, Kosove. 

Submitted January 6, 2015. Published January 2, 2008. 


1 





2 


M.A. ALQUDAH 


_ area(p 1 ,p 2 ,p) j nc [i ca te tb e barycentric ’’area” coordinates of the point p with 
respect to the triangle. 

Although there are three coordinates, there are only two degrees of freedom, since 
u + v + w = 1. Thus every point is uniquely defined by any two of the barycentric 
coordinates. That is, the triangular domain defined as 

T = {(it, v,w) : u, v, w > 0, u + v + w = 1}. 

Definition 1.1. The n + 1 Bernstein polynomials of degree n are defined by 

for i = 0,1 ,..., n, ( 1 . 1 ) 

where (™) is the binomial coefficients. Forf = (i,j,k) denote triples of non-negative 
integers such that |£| = i+j+k, then the generalized Bernstein polynomials of degree 
n are defined by the formula 

B™(u,v,w) = |C| = n, where ^ 

The generalized Bernstein polynomials have a number of useful analytical and 
elegant geometric properties [T[. In addition, the Bernstein basis polynomials of 
degree n form a basis for the space Tl n of polynomials of degree at most n. Thus, with 
the revolt of computer graphics, Bernstein polynomials on [0,1] became important 
in the form of Bezier curves, and the polynomials determined in the Bernstein 
(Bezier) basis enjoy considerable popularity in computer-aided Geometric design 
(CAGD) applications. 

Degree elevation is a common situation in these applications, where polynomials 
given in the basis of degree n have to be represented in the basis of higher degree. 
For any polynomial p(u, v, w ) of degree n can be written using Bezier coefficients 
d£ in the Bernstein form 

p(u,v,w) = ^ d c B2(u,v,w). (1.2) 

|C|=n 

With the use of degree elevation algorithm for the Bernstein representation, 
the polynomial p(u,v,w) in ( 11 . 21 ) can be written (multiplying both sides by 1 = 
u T v T w) as 

p(u,v,w ) = ^2 B™ +1 (u,v,w). 

I£l=”+i 

The new coefficients defined in [HIT] as d j. 1 '* = + kdij^—i) 

where |d = n+ 1. Moreover, the next integration is one of the interesting analytical 
properties of the Bernstein polynomials B™(u,v,w); 

JJ T B Z(u, v ,w) dA = 

where A is the double the area of T and (™J 2 ) is the dimension of Bernstein 
polynomials over the triangle. This means that the Bernstein polynomials partition 
the unity with equal integrals over the domain; in other words, they are equally 
weighted as basis functions. 
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Definition 1.2. Let p(u,v,w) and q(u,v,w) be two bivariate polynomials over T, 
then we define their inner product over T by 



With the inner product defined, we say that the two polynomials p(u, v, w) and 
q(u,v,w) are orthogonal if (p,q) = 0. 

For m > 1, let £ m denote the space of polynomials of degree m that are orthog¬ 
onal to all polynomials of degree < m over a triangular domain T, i.e., 

= {P ^ H-m . p _L 

For an integrable function f(u, v, w) over T, consider the operator S n (f) defined 
in as 

S„(/) = (n + l)(n + 2) £ (/, B™) B™. 

ICI=« 

For n> m, X m ,n = (n+m+$}n- m )\ is an eigenvalue of S n , and £ m is the corre¬ 
sponding eigenspace. The following lemmas will be used in the proof of the main 
results. 


Lemma 1.3. |5|. Let p = X^|c|=n C (B™ £ £ m and let q = X)|c|=n^C-®c e w 
m < n. Then, 

/ v __ ( n Q 2 _ p 

(n + m + 2 )!(n - m)\ ^ CC C ' 

Is I — n 

Lemma 1.4. [5] ■ Let p = ^|<j|= n C ( B ( e Then we have 

p £ £ n <*=> ^2 c ( d C = 0 Vg = ^2 d C B c e Tln-i- (1-3) 

|C|=n |C|=n 

1.2. Factorials. We present some results concerning factorials, double factorials, 
and combinatorial identities. The double factorial of an integer n is given by 

(2 n — 1)!! = (2 n — 1)(2 n — 3)(2n — 5)... (3)(1) if n is odd 

(1.4) 

n!! = (n)(n — 2)(n — 4)... (4) (2) if n is even, 

where 0!! = (—1)!! = 1. From (11.41) . we have the following definition. 


Definition 1.5. For an integer n, the double factorial is defined as 

( 2 2 ( 2 -)! if n is even 
if n is odd 


(1.5) 


From the definition, we can derive the factorial of an integer minus half as 


"- 5 " 


n!(2n — 

(2 n)H ' 


( 1 . 6 ) 


Moreover, the following identity can be used for the main results simplifications. 


Lemma 1.6. For an integer n, we have the following combinatorial identity 


,v 2 

n — k 


n — 


1 

2 ^ 


2 n 


(1.7) 
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Proof. By expanding the left-hand side and using m with some simplifications, 
we have 

(n-\\(n-\\ = (n- 1 )! (n-|)! 

\n—k)\ k ) (n — k)\(k — |)! k\(n — k — |)! 

(fin — 1 )!! 2 k ( 2 n — 1 )!! 2 n ~ k 
~ 2 n (n - k)\ ( 2 k - 1 )!! 2 "fc! ( 2 n - 2 fc - 1 )!! 

1 ( 2 n — 1 )!! ( 2 tt. — 1 )!! 

“ 2 n (n - k)lk\ (2k - 1 )!! ( 2 (n - k) - 1 )!!' 

Using the fact (2n)\ = (2 n — l)!!2"n! we get the desired identity. □ 


1.3. Univariate Chebyshev-I Polynomials. The Chebyshev-I polynomials T n (x) 
of degree n are the orthogonal, except for a constant factor, on [— 1 , 1 ] with respect 
to the weight function W(x ) = —s=y. 

The following two lemmas will be used in the construction of the bivariate 
Chebysliev-1 weighted orthogonal polynomials and the proof of the main results. 
The Pochhammer symbol is more appropriate, but the combinatorial notation gives 
more compact and readable formulas, these have also been used by Szego El- 


Lemma 1.7. [9] The Chebyshev-I polynomials T r (x) have the Bernstein represen¬ 
tation: 


T r (x) = 


r!(2r — 1)!! /2r 


2 2 r (2r)!! \r 


E 

2=0 


(-l) r 


i d(r-i)! 
(2z)!(2r — 2*)! 


Bf(x), r = 0,1, 


( 1 . 8 ) 


Lemma 1.8. [9] The Chebyshev-I polynomials Tq(x), ... ,T n (x) of degree less than 
or equal to n can be expressed in the Bernstein basis of fixed degree n by the following 
formula 

n 

T r (x)=Y J M™ r B?{x), r = 0,l,...,« 

2—0 


where 


M? r 



_1 min(2,r) 

E 

fc=max(0,i+r—n) 




(1.9) 


For simplicity and without lost of generality, we take x £ [0,1] for both Bernstein 
and Chebyshev-I polynomials. 


2. Chebyshev-I Weighted Orthogonal Polynomials 

In this section, a simple closed-form representation of degree-ordered system of 
orthogonal polynomials I7n^ r {u,v,w) is constructed on a triangular domain T, by 
generalization of the construction in [S]. These polynomials will be given in Bern¬ 
stein form, since Bernstein polynomials are stable [ 6 ]. The basic idea in this con¬ 
struction is to make I7 n , r (u,v,w) coincide with the univeriate Chebyshev-I polyno¬ 
mial along one edge of T, and to make its variation along each chord parallel to that 
edge a scaled version of this Chebyshev polynomial. The variation of ^ n ^ r {u,v,w) 
with w can then be arranged so as to ensure its orthogonality on T with every 
polynomial of degree < n, and with other basis polynomials 3? n ,s{'u>, v, w) of degree 
n for r ^ s. 









BIVARIATE CHEBYSHEV ON SIMPLICIAL 


5 


Lemma 2.1. [5]. For r = 0,..., n define the polynomials 

=zVi) j '( n+ * +1 W>)> (2- 1 ) 

j =0 \ 3 / 

then for i = 0,..., n — r — 1, Q n ,r{w ) is orthogonal to (1 — w ) 2r+l+1 on [0,1], and 
hence 

Qn,r(w)p{w)( 1 — -io) 2 r+ 1 d«; = 0 

/or every polynomial p{w) of degree less than or equal n — r — 1. 


Now, for r = 0,1,..., n and n = 0,1,2,... we define the bivariate polynomials 

+ r + 1 / 


£T nr {u,v,w) ='^2,c(i)B r i (u,v) 

2—0 j =0 


3 


Bf- r (w,u + v), (2.2) 


where 13/ (tt, u) is the Bernstein polynomials defined in (ED and 

i = 0,..., r. 


(2r\ /2r\ 

c(i) = (_l)r-*VrA2^ 


2 2 T) 


(2.3) 


To show that the bivariate polynomials 3n jr (u,v,w), r = 0,1,...,n; n > 0 
with respect to the weight function W(it, v, w), form an orthogonal system over the 
triangular domain T, we prove 3’ ntr (u,v,w ) £ £„, r = 0,1,..., n; n > 1, and for 
f -L '^n,s- 

Let ^ 0,0 = I, the polynomials 3n^ r (u,v,w) for r = 0,...,n; n > 0 form a 
degree-ordered orthogonal sequence over T. 

The polynomials in ( 12 . 21 ) can be written using the univariate Chebyshev-I poly¬ 
nomials form as: 


,r(«, W) = £ c(i)S^(l - < £(-l) J ( n + ^ + ^ Sr r ( W > 1 - «>)• 
»=o v ”r 1 j=o \ 3 / 


Using Lemma H771 and B\{u,v) = (u + v) r B^ (jff—), we get 


(2.4) 


2 2 t '(Y 1 ) 2 u 

3n !r (u,v,w) = , r r (--)(1 - w) r Q ni r(w), r = 0,... ,n, (2.5) 

( 2 r)! 1 — w 

where T r (t) is the univariate Chebyshev-I polynomial of degree r and <3 n , r (to) is 
defined in equation (12.11) . For simplicity and without loss of generality we rewrite 
(12.51) as 

u 

3n,r{u,V,w) = T r (- -) (1 - w) r Q n r (to) , T = 0, . . . , Tl. (2.6) 

1 — to 

Now, we show that the polynomials 3n t r(u,v,w), r = 0 ,n, are orthogonal 
to all polynomials of degree < n over the triangular domain T. 

For each s = 0,..., m and m = 0,..., n — 1 we define the bivariate polynomials 

g s , m (u, v, to) = T s ( U )(1 - w) m w n - m -\ m = 0,..., n - 1, s = 0,..., m. (2.7) 
1 — w 

The span of g s ,m{ tt, v, w) includes the set of Bernstein polynomials B™(u, v)w n ~ m ~ 1 . 
0,..., to; m = 0,..., n — 1, which span II n _i. 
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So, it is sufficient to show that for each s = 0,..., m; m = 0,... , n — 1, 

I := JJ X l , r (u,v,w)g Sjrn (u,v,w)W(u,v,w)dA = 0. ( 2 . 8 ) 

T 

The integral / can be simplified as 

1 1 — W 

I = a[ [ T r (-^—)Q nr (w)T s l—^—)w n - m - 1 u-iv-i(l-wV +r+m dudw. 

J J 1 — w 1—w 

0 0 

(2-9) 

By making the substitution u = t(l — w), we have 

l l 

I = A J T r (t)T a (t)t~$ (1 - t )~’ dt j QnA w )(. 1 - wy +r+rn w n - m - l dw. 

0 0 

If to < r then we have s < r, and the first integral is zero by the orthogonality 
property of the Chebyshev-I polynomials. If r < m < n — 1, the second integral 
equals zero by Lemma 12.II Thus we have the following theorem. 

Theorem 2.2. For each r = 0, l,...,n; n > 1, Ai,r{u,v, w) £ £„ with respect to 
the weight function W (u, v, w ) = 1 — w ) 7 such that 7 > 1 . 

Note that taking W (u, v, w ) = u~?v~A 1 — w ) T enables us to separate the inte¬ 
grand, and taking the constrain 7 > 1 enables us to use Lemma 12.11 

Now we need to show Ai,r{u, v, w ) is orthogonal to each polynomial (u, v, w) 
where r 7 ^ s. 


I = 


JJ &n,r{u, V , w)& niS (u, V , w)W(u, V , w)dA 


IT 
1 1—w 


= A 


T r 


0 0 


T„ 


1 — w J VI — w 


(1 - w) r+S Qn,r{w)Qn :S (w)W(u, V , w)dudw. 


By making the substitution u = f(l — w), we have 


1 1 

I = A j T r {t)T s {t)t~^ (1 - t)~* dt J Q nr (w)Q n ^ s (w)(l — w)' r+r+s dw. 

0 0 

The previous integral equals zero by orthogonality property of the Chebyshev-I 
polynomials, thus we have the following theorem. 

Theorem 2.3. For r ^ s, F? n Au,v,w) _L t% rljS (u,v,w) with respect to the weight 
function W (u, v, w ) = u~A~^ (1 — w) 7 such that 7 > —1 . 

Therefore, the bivariate polynomials Ai,r{u, v , re), r = 0,1,..., n; n > 0 form an 
orthogonal system over the triangular domain T with respect to the weight function 
W(u, v, w), 7 > —1. 
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3. Applications 


The Bernstein-Bezier form of curves and surfaces have some interesting geometric 
properties 0E1, which is very important in the numerical computations. So, the 
orthogonal polynomials ^ n ^{u,v,w), r = 0,1 ,..., n; n > 0 can be written in the 
following Bernstein-Bezier form: 

%,, r (u,v,w)= Y, a n c ’ r B?(u,v,w). (3.1) 

|C|=n 

We are interested in finding a closed form of the Bernstein coefficients a"’ r and 
derive a recursion relation that allow us to compute the coefficients efficiently. 
From equation (12.2D . ^ n ^ r (u,v,w) has degree < n — r in the variable w, so 

a ijk = 0 f° r k > n-r. (3.2) 

For 0 < k < n — r, the remaining coefficients are determined by equating m and 
EH as follows 


E <f t B~ k (u,v,w) = (- l)*r + ( +1 

i-\-j=n—k x 

Comparing powers of w on both sides, we have 


B k r (w,u + v)Y c ( i ) B i( u > v )- 


2=0 


n—k 


\ ' n ,r n ■ 

2 L %k „■ I1. 


i\j\k\ 


= (-!)* 


n + r + 1\ (n — r 


(u + v) n - r - k Y<i) B Uu,v). 


i =0 J x / \ / 2=0 

The left hand side of the last equation can be written in the form 


n—k 

E< 

2=0 

Now, we get 

n—k 


n,r n\(n-k)\ 

i\(n — k — i)\k\(n — k)\ 


n—k / \ 

2=0 X / 




n + r + 1\ /n — r 


(u+v) n - r - k Y^) B Hu,v). 


2=0 N 7 X / \ / 2=0 
With some binomial simplifications and using Lemma 11.81 we get 

n—k 


= (-!>*(■ 


n + r + 1\ (n — r 


Y M ?~ kB r k (u,v), 


2=0 


(3-3) 

where M™ r k are the coefficients resulting from writing Chebyshev-I polynomial of 
degree r in the Bernstein basis of degree n — k, as defined by expression (11.91) . Thus, 
the required Bernstein-Bezier coefficients given by the following theorem. 


Theorem 3.1. The Bernstein coefficients a"’ r of equation (EH are given explicitly 
by 


where Mf r k 


n,r _ 

a ijk ~ 


(-l) k 

0 


( + k + )( k ) A,rri-k 

(!) Mi ’ r 


given in (fL9l) . 


0 < k < n — r 
k > n — r 


(3.4) 
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To derive a recurrence relation for the coefficients a”’/ of J?n,r{u,v,w), consider 
the generalized Bernstein polynomial of degree n — 1, 


B ijk ( u , w ) = 


(i + 1) 


B i+i,j,k( u ’ v ’ w ) + 


(i + i) 


/i 


»j+i 

n—1 


,fc( 


x (& + 1) 

w, u, w)+- - -B 


i,j,k -\-1 


(u,v,w). 


From the construction of £? ntr (u,v,w), we have (BV\ k (u,v,w), 3? n ^ r {u,v,w)) = 0, 
where i + j + k = n— 1. Using Lemma. 1 1.41 


(3.5) 


By Theorem l3.il we have a"^_ i 0 = M" r for * = 0,1,..., n. Thus, we can use m 
to generate aU r fc recursively on k. 
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